New algorithm of the high-temperature expansion for the Ising model 
in three dimensions 



H. Arisuc^* and T. Fujiwara'' ^ 

^Osaka Prefectural College of Technology, Saiwai-cho 26-12, Neyagawa, Osaka 572-8572, Japan 

'^Faculty of General Studies, Kitasato University, Kitasato 1-15-1, Sagamihara, Kanagawa 228-8555, 
Japan 

New algorithm of the finite lattice method is presented to generate the high-temperature expansion series of 
the Ising model. It enables us to obtain much longer series in three dimensions when compared not only to the 
previous algorithm of the finite lattice method but also to the standard graphical method. It is applied to extend 
the high-temperature series of the simple cubic Ising model from fj'^^ to /3*® for the free energy and from /3^^ to 
for the magnetic susceptibility. 



1. INTRODUCTION 

The finite lattice method ||l|,|,| IS a power- 
ful tool to generate the exact high- and low- 
temperature series and other exact expansion se- 
ries for the spin models in the infinite volume 
limit. In the graphical method, one has to list 
up all the relevant diagrams and count the num- 
ber they appear. In the finite lattice method we 
can skip these jobs and reduce the main task to 
the calculation of the partition function for the 
finite size lattices, which can be done efficiently 
using the site-by-site integration without the 
graphical technique. 

It has been extremely effective primarily in two 
dimensions, but unfortunately it has worked out 
the series only in some limited cases in three di- 
mensions. This is because in two dimensions the 
CPU time and the computer memory needed to 
obtain the series to order N increase exponen- 
tially with N, while they grow up exponentially 
with iV^ in three dimensions. Here we present 
a new algorithm of the finite lattice method in 
which the CPU time and the computer mem- 
ory increase exponentially with NlogN. It en- 
ables us to generate much longer series of the 
high-temperature expansion for the Ising model 
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in three dimensions when compared not only to 
the previous algorithm of the finite lattice method 
but also to the graphical method. 

Although our main target is in three dimen- 
sions, the new algorithm of the finite lattice 
method applies in arbitrary dimensions, so we de- 
scribe it in two dimensions for convenience in the 
next two sections. 

2. FINITE LATTICE METHOD 

In the finite lattice method to generate the 
high-temperature series for the free energy in 
two dimensions we calculate the partition func- 
tion Z{lx X ly) for the finite size lattices with 
2{l.j; + ly) < N and define recursively |2| 

Hlx X ly) = iog[Z{l^ X ly)] 

E -^(Cxg. (1) 

Here we use the notation for the lattice size such 
that the 1x1 lattice means the unit square com- 
posed of 2 X 2 sites. The Boltzmann factor for 
each bond is expressed as 

exp iPskSk') = cosh(/3) (1 -I- tskSk') , (2) 

with (3 = J/ksT and t = tanh(/3). We define 
the bond configuration as the set of bonds to 
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which the factor tskSk' in is assigned while 
the factor 1 is assigned to the other bonds. Non- 
vanishing contribution to the partition function 
comes only from the bond configuration in which 
the bonds form one or more closed loops. Each 
of the closed loops is a polymer in the standard 
cluster expansion!^. Then the Taylor expansion 
of (p{lx X ly) with respect to t includes the con- 
tribution from all the clusters of polymers in the 
standard cluster expansion that can be embedded 
into the lattice of l^; x ly but cannot be embedded 
into any of its rectangular sub-lattices l'^ x ly and 
it starts from the term with n = 2{lx + ly), 
which comes from the cluster of a single polymer 
(one closed loop of bonds) that have two intersec- 
tions with any line perpendicular to the lattice 
bonds. The expansion series of the free energy 
density in the infinite volume limit is given by 



3. NEW ALGORITHM 



(3) 



In the standard algorithm of the finite lattice 
method the full partition function for the finite 
size lattice is calculated with all the bond con- 
figurations taken into account. In order to ob- 
tain the series to a given order, however, it is 
enough to consider only a restricted number of 
the bond configurations. Let us consider the 
anisotropic model of the simple cubic Ising model 
with (3i — Ji/kgT and ti ~ tanh(/3i) (i — x,y). 
To obtain the series for (l){lx x ly) to order Ny = 
2ly + ANy in ty we introduce [0 in the new algo- 
rithm (j){lx X ly,ANy) defined recursively by 



(j){l^ X ly,ANy) ^ log [Z{1^ X ly,ANy)] 
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cj,{l',xl'y,ANy) 



(4) 



Here the partition function Z{la: x ly, ANy) is cal- 
culated only with the bond configurations taken 
into account that have orders Uyi in ty for the i-th 
layer perpendicular to the y-direction satisfying 



E 



max(nyi, 2) 



< 2L 



AN,, 
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We neglect every bond configuration for the 
partition function that has ^ -^^^^ max(nj^i, 2) > 
2ly + ANy among the configurations that have 

X^iLi "-j/j '^ly + ANy. It is easy to prove that 
any of the neglected configuration does not con- 
tribute to X ly) in the order lower than or 
equal to Ny — 2ly + ANy. For such a configura- 
tion at least one of the n^i's should be zero, so 
they are disconnected configuration(composed of 
more than one polymer) or they can be embed- 
ded into a rectangular sub-lattice of l'^ x I'y with 
I'y < ly and in either case they do not contribute 
to 0(^2, X ^t/) in the order lower than or equal to 
Ny = 2ly + ANy. Thcy contribute to (t>{lx x ly) 



only in higher order than Ny 
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ANy by 



constituting the connected cluster of polymers to- 
gether with the polymers coming from other con- 
figurations that have n' y^ > 2 for the layer i with 
Uyi = . Examples of the bond configurations 
are shown in Figure for Ix = A, ly = 5. The 
example (a) has {uyi} — {2, 2, 4, 2, 2} and should 
be taken into account for ANy = 2, while the ex- 
ample (b) has {uyi} — {0,4,4,0,2} and should 
be neglected for the same ANy = 2, in spite of 
the fact that the total order of (b) in ty is smaller 
than 2ly + ANy. 



(a) 



(b) 



Figure I. Examples of the bond configurations. 



The contribution of the bond configuration 
with {fiyi} to the partition function of the finite 
size lattice can be calculated by the transfer ma- 
trix formalism as 



z{M) 



(6) 



Here Vj-j-_^-^^ is the transfer matrix element with 
incoming Uyi spins and outgoing Uy spins. 
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The summations over the spin locations ji , j2 , • ■ ■ 
of the ?T.j,i, nj,2, • ■ • spins, respectively, are as- 
sumed in the right hand side of @. 

In three dimensions, which is our main target, 
the transfer matrix element itself can be regarded 
as the partition function in two dimensions and 
can be calculated efficiently by the site-by-site 
construction ^,||. In the new algorithm most of 
the CPU time should be used for the calculation 
of these transfer matrix elements, and the total 
CPU time to generate the series to order N can be 
estimated to increase exponentially with the lead- 
ing term of the exponent proportional to N log N. 

4. RESULT 

The new algorithm is applied to generate the 
high-temperature series of the simple cubic Ising 
model for the free energy Q to P"^^ and for the 
magnetic susceptibility § to The obtained 

series agree with the previous series to for the 
free energy which was obtained by the previ- 
ous algorithm of the finite lattice method, and 
to /3^^ for the susceptibility [p^ , which was ob- 
tained by the graphical method. It should be 
commented that the previous algorithm of the fi- 
nite lattice method can generate the susceptibility 
series only to P^^. 

Preliminary analysis of the free energy series 
using the inhomogeneous differential approxima- 
tion and the ratio method biased by the value 
of the critical point f3c gives the estimation of 
the critical exponent for the specific heat as 
a = 0.104(1) and a = 0.108(1) correspond- 
ing to the result of f3c = 0.22165459(10) and 
Pc = 0.2216595(15), respectively, of the recent two 
Monte Carlo simulations pT| , p^ . Unfortunately 
the free energy series can give the estimate of the 
critical point itself only in poor precision. The 
analysis of the susceptibility series by the ratio 
method gives the estimate of the critical point 
as Pc = 0.2216550(5) and the critical exponent 
as 7 = 1.2370(2). This estimation for the criti- 
cal exponent do not use the value of the critical 
point as the input. It is consistent with the recent 
estimate 7 = 1.2371(4) obtained from the high- 
temperature series of the generalized Ising model 



5. DISCUSSION 

We have presented the new algorithm of the fi- 
nite lattice method for the high-temperature ex- 
pansion of the Ising model. It has been applied to 
the high-temperature expansion of the free energy 
and the magnetic susceptibility for the simple cu- 
bic Ising model. It can be applied to the high- 
tcmperature expansion of other quantities such as 
the correlation length and it can also be applied 
to the models with continuous spin variables such 
as the XY model in three dimensions. We note 
that the dimensionality of the lattice is not re- 
stricted to three as can seen by the fact that the 
description of the new algorithm in section 3 was 
given in two dimensions. Furthermore the basic 
idea of the new algorithm can be used in the low- 
temperature expansion for the spin models. We 
can expect that the new algorithm will enable us 
to generate the series for these models that are 
much longer than the presently available series. 
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